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Abstract. Accretion-powered X-ray pulsars are among the most luminous X-ray sources in the 
Galaxy. However, despite decades of theoretical and observational work since their discovery, no 
satisfactory model for the formation of the observed X-ray spectra has emerged. In this paper, we 
report on a self-consistent calculation of the spectrum emerging from a pulsar accretion column that 
includes an explicit treatment of the bulk and thermal Comptonization occurring in the radiation- 
dominated shocks that form in the accretion flows. Using a rigorous eigenfunction expansion 
method, we obtain a closed-form expression for the Green's function describing the upscattering 
of monochromatic radiation injected into the column. The Green's function is convolved with 
bremsstrahlung, cyclotron, and blackbody source terms to calculate the emergent photon spectrum. 
We show that energization of photons in the shock naturally produces an X-ray spectrum with a 
relatively flat continuum and a high-energy exponential cutoff. Finally, we demonstrate that our 
model yields good agreement with the spectra of the bright pulsar Her X-l and the low luminosity 
X Per. 

Keywords: Neutron stars, Accretion, X-Rays, Comptonization, Radiation-dominated shocks 
PACS: 95.30.Jx, 97.60.Gb, 97.80.Jp, 98.35.Mp 

INTRODUCTION 

Accreting X-ray pulsars were discovered in the late 1960s by rocket flights from White 
Sands in New Mexico yj,|2|]. Since that time the number of known X-ray pulsars of all 
types has grown to over 100. The prevailing model for accreting X-ray pulsars consists of 
a neutron star orbiting a normal stellar companion that is losing matter via either Roche 
lobe overflow or a stellar wind. Accreting X-ray pulsar spectra are characterized by a 
power-law dependence on energy in the range above ~ 3 keV with a quasi-exponential 
cutoff at higher energies, typically near 20-40 keV. For example, see the X-ray spectrum 
of 4U 11626-67 reported by Orlandini y|], or the multiple source spectra reported by 
Coburn et al. [4]. Cyclotron features have also been observed in many sources |4J]. 

Comptonization is known to play an important role in the X-ray pulsar spectral forma- 
tion process. This is a general result, based on the overall spectral shape (a power-law) 
and the fact that Comptonization tends to result in a power-law dependence with energy 
in many accreting compact object sources. The usual method of characterizing accreting 
X-ray pulsar spectra is to fit the broad-band X-ray spectra using a number of ad-hoc 
functions such as power-laws, Gaussian emission line features, Gaussian or Lorentzian 
absorption features, and various types of quasi-exponential high-energy cutoffs. How- 
ever, most of these functional forms have little physical motivation beyond the fact that 
they "look like the spectra" and they are easy to incorporate into spectral fitting pro- 



grams (e.g., XSPEC). The real physical parameters of the source (e.g., accretion rate, 
accretion region size, shock height, plasma temperature, etc.) are not directly connected 
with the ad hoc parameters and are in fact much harder to determine. 

Three principal models have been put forward to explain the energy spectra of accret- 
ing X-ray pulsars: the gas-mediated collisionless shock model [HI]; the Coulomb colli- 
sional stopping model Bo, 13, ISD ; and the radiation-dominated flow model §, EU. The 
collisionless shock model and the Coulomb collisional model can only be applied in the 
case of low-luminosity accretion onto neutron stars. In such models, the effects of radi- 
ation pressure are assumed to be small, and the flow impinges directly onto the neutron 
star surface (in the absence of other effects; see below). These models tend to produce 
pencil-beamed emission patterns because the radiation escapes primarily through the 
top of the neutron star atmosphere. In high-luminosity pulsars, the dynamical structure 
is expected to be dominated by the effects of radiation pressure, which decelerates the 
gas to rest at the stellar surface. The radiation-dominated inflow models of Davidson 
Q9U and Arons et al. jl(X\ attempt to describe the flow dynamics across a broad range of 
X-ray pulsar luminosities, all the way up to the Eddington limit. 

Neither the collisionless shock model nor the Coulomb collisional model have demon- 
strated good agreement with actual X-ray pulsar spectra. Meszaros and Nagel [11] com- 
pared the spectrum of Her X-l with results obtained using the Coulomb collisional 
model. The observed Her X-l spectrum was not well fit by the calculated spectra. 
Furthermore, the applicability of the Coulomb collisional model to Her X-l is ques- 
tionable because the source luminosity is believed to be close to the Eddington limit 
(L x ~ 2.2 x 10 37 ergs s _1 at 5 kpc), implying that radiation pressure is important 



[|12L 11311 . For further discussion of Coulomb collisional stopping models see Harding 
and Lai [[bill and references therein. 

Insight into the photon transport in the accretion column can be gained by estimating 
the optical depth to electron scattering in the magnetic field direction, Tii, and the 
scattering optical depth across the accretion column, Tj_. Assuming Thomson scattering 
and a free-fall velocity profile, we can express these quantities in cgs units using [see 
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where L x is the accretion luminosity, M ns and R ns are the neutron star mass and radius, 
respectively, and tq is the radius of the accretion region on the neutron star surface. Based 
on the modeling described below, we find that Tii ~ 6.5 x 10~ 2 and Tj_ ~ 4.2 x 10~ 3 for 
X Per. Conversely, in the Her X-l case we obtain Tji ~ 2.1 x 10 4 and t_l ~ 1.4 x 10 2 . 
Hence photons will scatter many more times while escaping from the accretion flow 
in Her X-l than in X Per. This is a general result that applies to all high-luminosity 
X-ray pulsars such as Her X-l because of the larger accretion rates compared with the 
low-luminosity, steep spectrum sources such as X Per. 



A RADIATION-DOMINATED SHOCK MODEL FOR X-RAY 

PULSAR SPECTRA 



Adopting a cylindrical, plane-parallel geometry for the accretion column with the mag- 
netic field in the z-direction, the Green's function f G (zo,Z, £o,e) satisfies the modified 
[fl6| steady-state transport equation [see, e.g.. fl7lfl8|l 
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where z is the altitude above the stellar surface, v < is the inflow velocity, Nq is the 
rate of injection of seed photons with energy £o at location zo, t esc represents the mean 
time photons spend in the plasma before diffusing through the walls of the column, a;, is 
the electron scattering cross section for photons propagating parallel to the magnetic 
field, <7 is the angle-averaged cross section, and T e , n e , and m e denote the electron 
temperature, number density, and mass, respectively. The mean escape time is computed 
using f esc = ro/w±, where = c/t± is the diffusion velocity perpendicular to the z- 
axis, t_l = n e o±ro is the electron scattering optical thickness across the column, and 
o± denotes the electron scattering cross section for photons propagating perpendicular 
to the magnetic field. The solution for the Green's function / G (zo,z, £o 5 e) is obtained 
by deriving eigenvalues and associated eigenf unctions based on the set of spatial and 
energetic boundary conditions for the problem [seeO]. We define a source function 
Q such that e 2 Q(z,e)dedz gives the number of seed photons injected per unit time 
in the altitude range z to z + dz and energy range £ to e + de. Once we have the 
analytical solution for the Green's function, the particular solution corresponding to 
bremsstrahlung, cyclotron, or blackbody source distributions can be obtained via the 
integral convolution 

/fee) = r rtd^pA e 2 Q(zo,£o)deodz . (4) 
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In the models described here, the effects of bulk and thermal Comptonization are treated 
explicitly using the transport terms in equation ©. The importance of dynamical (bulk) 
versus thermal Comptonization depends on the parameter 5 = (aa || /3a)(m e c 2 /kT e ), 
where a ~ 0.3 — 0.5 describes the velocity variation as a function of the optical depth 
above the stellar surface. It can be shown that 8 is essentially the ratio of the "y- 
parameters" for bulk and thermal Comptonization. When 8 is of order unity, the two 
processes are comparable, and when 8 ^> 1, the bulk process dominates. Another im- 
portant spectral formation parameter is ^ = (nrom p c)/(M ^o^o±), where M denotes 

the accretion rate and m p is the proton mass. We find that t, is roughly equal to the ratio 
of the dynamical (accretion) timescale divided by the timescale for the photons to dif- 
fuse through the column walls. The condition ^ ~ 1 must be satisfied in order to ensure 
that radiation pressure decelerates the gas to rest at the stellar surface. 



RESULTS FOR X-RAY PULSAR SPECTRA 



Using our model we can compute the theoretical spectrum emitted from an X-ray pulsar 
accretion column due to Comptonized bremsstrahlung, cyclotron, and blackbody seed 
photons. The theoretical phase-averaged photon count rate spectrum, F £ (e), is given by 

F £ (e) = ^§, * B (e)= r^f( z ,e)dz, (5) 
47ix> / Jo t esc (z) 

where 4> e (e) represents the vertically-integrated escaping photon number spectrum, 
f(z, e) is computed using equation ©, and D is the distance to the source. As a check 
on our results for the spectra, we confirm that the number of photons escaping from 
the column per unit time is exactly equal to the number injected, as required by our 
steady- state scenario. 

In Figured] we plot the theoretical count-rate spectrum F E (e) evaluated using equa- 
tion © along with the deconvolved, phase-averaged BeppoSAX spectrum of Her X- 1 
reported by dal Fiume et al. Ill 911 . The theoretical parameters in this case are a = 0.40, 
% = 1.45, Oj_ = ct t , 8 = 1.8, B = 3.80 x 10 12 G, M = 1.11 x H^gs -1 , r = 44m, 
T e = 6.25 x 10 7 K, and 7^ = 5.68 x 10 7 K, where o T denotes the Thomson cross section 
and T t h is the thermal mound temperature. We assume a source distance of D = 5kpc. 
In Figure \T\ results are plotted for the total spectrum, as well as for the individual con- 
tributions to the observed flux due to the Comptonization of cyclotron, blackbody, and 
bremsstrahlung seed photons. The theoretical spectrum in FigureCQalso includes an iron 
emission line. 

The general shape of the theoretical spectrum agrees with the observations for Her X- 
1 quite well, including both the quasi-exponential cutoff energy and the power-law slope. 
In the case of Her X- 1, reprocessed (Comptonized) blackbody emission from the thermal 
mound makes a negligible contribution to the spectrum because the radius of the accre- 
tion column is relatively small. Due to the high temperature of the post-shock plasma, the 
reprocessed cyclotron emission is overwhelmed by reprocessed bremsstrahlung emis- 
sion, which dominates the observed spectrum. 

The second source, 4U 0352+30 (X Per), is a low-luminosity (L x ~ 10 34 ergs s _1 ) 
X-ray pulsar with a relatively steep spectrum. This pulsar was discussed by Becker and 
Wolff B, l21ll using a pure-bulk Comptonization model, and it therefore represents an 



interesting test for the thermal+bulk model developed by Becker and Wolff QJ/ZU and 
described in this paper. In Figure [2] we compare the theoretical spectrum with archival 
RXTE data taken in July 1998 and reported by Delgado-Marti et al. I122n . The theory 
parameters in this case are a = 0.51, | = 1.85, a± = o T , 5 = 10.9, B = 3.30 x 10 12 G, 
M = 3.23 x 10 13 gs-\ r = 430m, T e = 4.00 x 10 7 K, T th = 9.00 x 10 6 K, and D = 
0.35kpc. The distance and magnetic field values are from Negueruela ll23ll and Coburn 
et al. ||24|1 . respectively. In contrast to the case of Her X-l, our model shows that the 
spectrum of X Per is dominated by Comptonized blackbody emission, which is due 
to the order of magnitude increase in the radius of the accretion column. Reprocessed 
cyclotron and bremsstrahlung radiation make a negligible contribution to the observed 
spectrum for this source. The spectral results obtained here using the thermal+bulk 
Comptonization model are nearly identical to those obtained by Becker and Wolff 
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FIGURE 1. Theoretical spectrum of Her X-l based on our radiation-dominated, radiative shock model 
(eq. J5|), compared with the data reported by dal Fiume et al. Il9ll . 
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FIGURE 2. Theoretical spectrum of X Per based on our radiation-dominated, radiative shock model 
(eq. (5)), compared with the data reported by Delgado-Marti et al. 112211 . Various amounts of interstellar 
absorption have been included as indicated. 



[|21i|20D, which is consistent with the fact that the spectral formation in this source is 
dominated by bulk Comptonization, as indicated by the large value of 8. 



CONCLUSIONS 



We have developed a new analytical model describing the spectral formation process 
in accretion-powered X-ray pulsars. The model includes a rigorous treatment of both 
the bulk and thermal Comptonization occurring in the radiation-dominated, radiative 



shock. These two types of Comptonization influence different regions of the radiation 
distribution and can explain a wide range of accretion-powered X-ray pulsar spectra. We 
have shown that the theoretical spectra produced by our model in luminous sources such 
as Her X-l are dominated by Comptonized bremsstrahlung emission, not Comptonized 
cyclotron emission as was previously conjectured. On the other hand, we find that 
the spectra of low-luminosity sources such as X Per are dominated by Comptonized 
blackbody radiation. Our new Comptonization model provides greatly improved fits to 
the observed spectral data when compared with the gas-mediated collisionless shock or 
Coulomb collisional models. Furthermore, the X-ray spectra produced by our model can 
be evaluated using a single-pass, closed-form algorithm that does not require numerical 
iteration. Our model should therefore be suitable for incorporation into the XSPEC 
spectral modeling environment and we expect to complete that work in the near future. 
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